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ABSTRACT 

In a previous work, we reported on the discovery of supersonic magnetic upfiows on granular cells in data from the 
SuNRiSE/IMaX instrument. In the present work we investigate the physical origin of these events employing data of the 
same instrument but with higher spectral sampling. By means of the inversion of Stokes profiles we are able to recover 
the physical parameters (temperature, magnetic field, line-of-sight velocity, etc) present in the solar photosphere at 
the time of these events. The inversion is performed in a Monte-Carlo-like fashion, that is, repeating it many times 
with different initializations and retaining only the best result. We find that many of the events are characterized 
by a reversal in the polarity of the magnetic field along the vertical direction in the photosphere, accompanied by an 
enhancement in the temperature and by supersonic line-of-sight velocities. In about half of the studied events, large 
blue-shifted and red-shifted line-of-sight velocities coexist above/below each other. These features can be explained 
in terms of magnetic reconnection, where the energy stored in the magnetic field is released in the form of kinetic 
and thermal energy when magnetic field lines of opposite polarities coalesce. However, the agreement with magnetic 
reconnection is not perfect and therefore, other possible physical mechanisms might also play a role. 

Subject headings: Sun: surface magnetism - Sun: magnetic topology - Sun: photosphere - Sun: 
granulation - polarization 



Section 6 briefly addresses the choice of model for the in- 
version. Finally, Section 7 presents our main conclusions. 



I. INTRODUCTION 

In the quiet solar photosphere the energy stored in 
the magnetic field is comparable to the kinetic energy 
due to convective motions. This gives rise to a rich 
variety of phenomena that evolve at very short time 
and spatial scales. The IMaX instrument (Imaging 
Magnetograph eXperiment; Martinez Fillet et al. 2011a 
) onboard of stratospheric balloon Sunrise (Barthol et 
al. 2011; Solanki et al. 2012) has helped to uncover 
many of these phenomena (see e.g. Solanki et al. 2010), 
from resolving magnetic flux-tubes (Lagg et al. 2010), 
to finding vortex tubes (Steiner et al. 2010), vortex 
flows (Bonet et al. 2010), etc. One of those discoveries 
involves supersonic magnetic upfiows (Borrero et al. 
2010, 2012). These events are characterized by highly 
blue-shifted circular polarization signals, that appear at 
the center or edges of granular cells and last for about 
80 seconds. Similar events were subsequently found also 
in data from the SF instrument onboard the Hinode 
spacecraft (Martinez Fillet et al. 2011b). The fact 
that magnetic fields of opposite polarity connected by 
horizontal fields appear in the vicinity of these events 
in 70 % of the cases, led us to surmise that they are 
caused by magnetic reconnection. In this paper we 
study them in more detail using a different data set 
from the IMaX instrument. A comparison between 
the old and new data sets is provided in Section 2. 
In Section 3 we describe the criterion employed to 
select events. The analysis technique, namely, the 
inversion of the observed Stokes profiles to retrieve the 
physical conditions of the solar atmosphere, is detailed 
in Section 4. Section 5 presents our results while 

^ ,—— this work, however, we will use a different observing 

|borrero@kis.uni-freiburg.de, vmp@ll.iac.es, woIfgang@kis.uni-freiburg .ifi)^f^p}jiiC^f^p,<j%] j^^€^l^i£[cqc^^ j^^^ intensity 



2. INSTRUMENTS AND OBSERVATIONS 

The data employed in this work were recorded with 
the stratospheric balloon-borne observatory Sunrise 
(Solanki et al. 2010; Barthol et al. 2011). Sunrise was 
launched on June 8, 2009 from Kiruna (Sweden) and 
landed on June 13, 2009 on Somerset Island (Canada). 
During this time. Sunrise's 1-meter telescope took 
broad-band images in different spectral windows with 
the SUFI instrument (Gandorfer et al. 2011), and 
spectropolarimetric data of the solar photosphere with 
IMaX (Martinez Fillet et al. 2011a). An average 
fiight altitude of 35 km allowed Sunrise to avoid more 
than 95 % of the disturbances introduced by Earth's 
atmosphere. In addition, image motions due to wind 
during the flight were stabilized by the Correlation- 
Tracker and Wavefront Sensor (CWS; Berkefeld et al. 
2011). Owing to the aforementioned advantages, IMaX 
spectropolarimetric data yielded a spatial resolution of 
0.25" and a field-of-view of 50"x50". Further image 
reconstruction based on phase diversity calibration of 
the FSF of the optical system improved the resolution 
to 0.15"-0.18". 

In Borrero et al. (2010) we employed reconstructed 
IMaX data that included the four components of the 
Stokes vector (/, Qj U , V) measured at five wavelength 
positions across the Fe I 5250.217 A spectral line. In the 
following we will refer to this observing mode as V5-6. 
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/ and circular polarization V were measured in twelve 
(instead of five) wavelength positions). For reasons that 
will be explained later, we restricted ourselves to employ 
non-reconstructed data with a spatial resolution of only 
0.25". 

Figure 1 displays the region of the intensity spectra 
around the Fe I 5250.217 spectral line, as recorded by 
the Fourier Transform spectrometer in the quiet Sun 
(Wallace et al. 1998). Crosses in this figure show the 
five wavelength positions scanned by the V5-6 data 
used in Borrero et al. (2010). These were located at 
A = [-80.-40,40,80,227] mA from line-center. Filled 
circles illustrate the wavelength positions scanned in 
the L12-2 observing mode. Here, the wavelength range 
goes from —192.5 mA to +192.5 mA, in twelve positions 
equidistantly distributed in steps of 35 mA. 
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Fig. 1. — Comparison of the V5-6 (crosses) and L12-2 (filled 
circles) observing modes. The former mode records the four com- 
ponents of the Stokes vector, while the latter acquires Stokes / and 
Stokes V (see text for details). The solid- black line corresponds to 
the Fourier Transform Spectrometer data (FTS-atlas) from Kitt 
PeaJc observatory in Arizona. The effective Lande factors g^s of 
eax;h line are calculated under the LS approximation from the elec- 
tronic configurations given in Table 1. 

The lack of linear polarization profiles Q and U 
in the L12-2 observing mode makes the polarimetric 
calibration of the data slightly more difficult to imple- 
ment compared to V5-6. The instrument's calibration 
matrix had been theoretically calculated (Martinez 
Fillet et al. 2011a) and experimentally confirmed at the 
INTA (Instituto Nacional de Tecnicas Aeroespaciales) 
facilities in Spain (Martinez Fillet 2007; del Toro Inicsta 
& Martinez Fillet 2012) stich that linear polarization 
cross-talk could be minimized by tuning the voltages 
of the nematic liquid crystals to the appropiate re- 
tardances. Unfortunately, this calibration stage was 
not performed while mated to the Sunrise telescope 
(main and secondary mirrors) which introduced more 
cross-talk than anticipated. Therefore, we expected to 
find a non-negligible contribution from Q and U in the 
measured linear combinations of / and V. Fortunately, 
the events that we will focus on, the so-called supersonic 
magnetic upflows in granular cells, had already been 



analyzed using the V5-6 observing mode. From that 
data we learned that these events have a negligible 
amount of linear polarization (although patches of 
enhanced linear polarization usually appear within 2" 
of these events; see for instance Fig. 4 in Borrero et al. 
2010 and Figs. 3-7 in Borrero et al. 2012). This means 
that, even in our scenario, where there is a significant 
cross-talk from Q and U into Stokes V, the circular 
polarization is not affected much. The uncorrected 
cross-talk does however slightly increase the noise level 
in the circular polarization. 

IMaX used the L12-2 mode to observe the solar 
photosphere for about 60 minutes on June 10, 2009. 
The observations were recorded close to disk center, 
/X = cos 9 = 0.99, where 9 corresponds to the heliocen- 
tric angle. Due to pointing problems, the time series 
was interrupted several times and only about half of 
that time is usable. In total, there are 52 full scans of 
Stokes I and V in twelve wavelength positions. Each 
scan is recorded during a time interval of 33 seconds. 
The noise-level is estimated to be 10~^ in units of the 
average quiet Sun continuum intensity. 



3. SELECTION OF EVENTS 

In Borrero et al. (2010), where we employed V5-6, 
the supersonic magnetic upflows were detected in the 

circular polarization (Stokes V) at AA = +227 mA 
(see Fig. 1). At this wavelength any signal can be 
produced by cither a strong red-shift (downflow) from 
Fe 1 5250.217 A or by strong blue-shift (upflow) from 
Fe I 5250.653 A. In the aforementioned work, the signal 
was ascribed to the latter case because the Stokes I 
signal from Fe I 5250.217 A was blue-shifted. 

Bearing this information in mind, we look for a 
strategy to find highly blue-shifted Stokes V profiles in 
the L12-2 data. The first idea that comes to mind is 
to use the blue-most wavelength position in these data, 
AA = —192.5 mA, to find such upflows. Unfortunately, 
this scanning position is not far enough towards the blue 
to guarantee that the observed signal at this wavelength 
is only due to large velocities, as it could also be 
caused by a magnetic field that shifts the cr-component 
of Stokes V into this wavelength. Indeed, displaying 
V{X = Ao — 192.5 mA), reveals a pattern that closely 
resembles the network. Thus, this wavelength alone 
cannot be employed to uniquely identify large upflows. 

In order to disentangle network elements from possi- 
ble supersonic magnetic upflows, we propose a different 
strategy, in which we compare the circular polarization 
close to the spectral line core and the circular polariza- 
tion close to the blue continuum. Let us refer to these 
two quantities as Vc and Viine) respectively. They are 
defined as: 



Mine = 



.1=4 



i=7 



(1) 
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Vc = ^Y.\^{Xi)\, (2) 
^ i=i 

where the index i runs from the blue- most {i = 1) o the 
red-most (i = 12) scanning positions indicated by the 
filled circles in Fig. 1. We note that, in the definition 
of l^ino we are subtracting Stokes V in the red wing 
{i = 6, 7, 8) from Stokes V in the blue wing (i = 3, 4, 5). 
This is done in order to obtain the polarity of the 
magnetic field vector, and has the additional benefit of 
partially canceling the noise. 

In the top-left panel of Figure 2 we display Vune 
(normalized to the averaged quiet Sun continuum 
intensity Iqs) over a portion of the full field-of-view 
from one of our available 52 snapshots. The regions 
of enhanced Vunc correspond mostly to the network 
elements since the circular polarization close to the 
line-center is large. In the top-right panel of Figure 2 
we plot, for the same region, the absolute value of the 
quotient of Vc and Mine- In this panel, the network 
appears as those regions where ||T4/Viinc|| ~^ 0. Regions 
where HVc/MincH >> 1 denote Stokes V profiles that 
are highly blue-shifted. Our selection criterion will 
consider as supersonic magnetic upflows any pixel in the 
field-of-view where ||K;/T4ine|| > 4. In Figure 2, those 
regions are indicated by the white contours. Note that 
the selected events also coincide with the center/edges 
of gramilar cells where the linc-of-sight velocity is 
blue-shifted by about —2 km s^^ (sec black contours 
in the left-bottom and right-bottom panels in Figure 
2). This confirms that the selected events have the 
same properties as those studied in Borrero et al. (2010). 

The linc-of-sight velocity Vlos displayed in Figure 2 
is obtained by calculating the center-of-gravity of Stokes 
I and correcting for gravitational red-shift, convective 
blue-shift, and for the wavelength shift across the FOV 
due to the collimated configuration of the instrument 
(see Sect. 9.1 in Martinez Fillet et al. 2011a). We 
mention this in order to avoid confusion with the Vlos 
that will be used in the next sections of this paper, 
which will be inferred from the simultaneous fitting of 
Stokes / and V. 

It is also important to clarify that Figure 2 shows a 
somewhat exceptional situation, in which three events 
occur on a small portion of the full field-of-view. This 
case has been selected to highlight the properties of 
the selected events, but by no means corresponds to 
the typical case seen in the observations. In fact, from 
the 52 available snapshots, the aforementioned selection 
criteria selects 857 pixels, belonging to 122 events. This 
results in an average of 2.3 events in each 50"x50" 
snapshot. Unfortunately, the pointing problems de- 
scribed in Sect. 2 prevented us from having a continuous 
time-series and therefore we cannot track each event in 
time. Thus, some of those 122 events might correspond 
to the same one but at different times in their evolu- 
tion. Consequently, we cannot compare these numbers 
with the occurrence rates obtained in Borrero et al. 2010. 



Once we have the selected pixels that correspond to 
the possible supersonic magnetic upflows in granules, we 
now proceed to extract physical information from their 
corresponding Stokes / and V profiles. This is done by 
means of the inversion of the radiative transfer equation 
employing the SIR (Stokes Inversion based on Response 
functions) code (Ruiz Cobo & del Toro Iniesta 1992). 
Starting with an initial model of the solar photosphere, 
SIR solves the radiative transfer equation to obtain the 
theoretical Stokes vector that arises from such model. 
The observed Stokes vector is then compared to the 
theoretical one through a merit-function. Via a 
Levenberg-Marquardt method, the original model is 
then iteratively modified until reaches a minimum. 
At each iteration step, the Levenberg-Marquardt al- 
gorithm provides the perturbations in the physical 
parameters at several optical-depth positions called 
nodes, that are needed to produce a better fit to 
the observed Stokes profiles. Each node represents a 
free parameter in the inversion. The resulting model 
(from the x^-minimization) can be then considered to 
represent the physical conditions present in the solar 
photosphere. For some recent reviews on this subject we 
refer the reader to del Toro Iniesta (2002), Bellot Rubio 
(2006) and Ruiz Cobo (2007). Since photon noise affects 
the results of the inversion rather negatively (see e.g. 
Borrero & Kobel 2011, 2012) and owing to the fact that 
image reconstruction techniques slightly increase the 
noise in the observations, we considered for the inversion 
the unreconstructed data with a spatial resolution of 
0.25 arcsec (see Sect. 2). 

Our inversions have been carried out with a 1- 
component model, in which the photosphere is consid- 
ered to be laterally homogeneous within each selected 
pixel. Thus, we only need to consider the vertical vari- 
ations of the physical parameters in the photosphere. 
These variations are often treated in terms of the dimen- 
sionless optical-depth at a reference wavelength of 5000 
A, T5, instead of the geometrical height z. In general, 
the physical parameters relevant for the formation 
of spectral lines are: temperature T{t^), line-of-sight 
velocity Vlos(''"5)! the three components of the 

magnetic field vector: B{t5) (modulus of the magnetic 
field vector), 7(r5) (inclination of the magnetic field 
vector with respect to the observer's line-of-sight) and 
4i{t^) (azimuth of the magnetic field vector in the plane 
perpendicular to the observer's line-of-sight). Other 
quantities, such as the gas Pgir^) and electron Peir^) 
pressure, as well as the density /^(ts), are derived from 
r(T5), the condition of hydrostatic equilibrium, and the 
equation of ideal gases with a variable mean molecular 
weight. In our invcrions we allow for the following free 
parameters (nodes): three for T(t5), one for B{tc,) (con- 
stant value with height), two for 7(t5), five for Vhosir^), 
and finally one for the micro-turbulent velocity 14iic(T5) 
(also constant with height). This adds up to a total of 
12 free parameters. The full stratification of the physical 
parameters with T5 is obtained via interpolation across 
the values at the nodes. In section 6 we give more 
details about our choice of model and free parameters, 
as well as discussing its implications. 



4. INVERSION OF STOKES PROFILES 
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Fig. 2. — Top-left: close-up (28" X 25") of a snapshot displaying the total circular polarization around the center of the spectral line Viine 
(Eq. 1) and normalized to the average continuum intensity over the quiet Sun Iqs- Top-right: same FOV as before but showing ||Vc/yiinclli 
where Vc corresponds to the circular polarization on the blue- wing of the spectral line (Eq. 2). Bottom-left: same FOV as before but showing 
the continuum intensity Ic normalized to the average continuum intensity over the quiet Sun /qs. Bottom-right: line-of sight velocity Vlos 
derived from the center-of-gravity of Stokes V. In all panels the black and white contours enclose the regions where ||Vc/Viincll > 4. These 
patches contain pixels that are selected for our study (see text for details). While the upper panels are unreconstructed, the bottom ones 
have been subject, for visualization purposes, to image restoration (see Sect. 2). 
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TABLE 1 

Atomic parameters of the spectral lines included in the 
inversion. 



200 



Specie 




Xlow 

[eV] 


log(9/) 


Elec.conf 


a 2 


a 2 




Co I 
Fe I 
Fe I 


5250.008 
5250.217 
5250.653 


4.175 
0.121 
2.198 


-0.114 
-4.938 
-2.198 


4/^ 4 rr 
'-'5/2 — -"7/2 

"■Do - ^Di 

_ 


n/a 
207 
344 


n/a 
0.253 
0.268 


0.785 
3.0 
1.5 



^ Aq represents the central wavelength position of the spectral line 
on the Sun. 

2 (T and a represent the atomic transition's cross-section (in units 
of Bohr's radius squared Oq) and velocity parameter, respectively, 
for collisions with neutral atoms under the ABO theory (Asntee 
& O'Mara 1995; Barklem et al. 1998). Collisional data for Co I 
is not available, and therefore we rely on Unsold's theory (Unsold 
1955) to calculate the collisional broadening of the spectral line. 

inverted a total of 100 times, where the initial values of 
the physical parameters are randomly chosen each time. 
From all those 100 independent inversions we retain only 
the one with the smallest value of x^- In the inversion 
we include the effects of the three spectral lines present 
in Figure 1. This is done because it is not possible to 
rule out the possibility of the Co I 5250.008 A and/or 
Fe 1 5250.653 A spectral lines entering the wavelength 
regions scanned by the L12-2 observing-mode (see filled 
circles in Figure 1), specially when considering events 
that involve large line-of-sight velocitie^. The atomic 
parameters for these three spectral lines are given in 
Table 1. 
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Fig. 3. — Observed Stokes / (upper panel) and Stokes V (bottom 
panel) profiles from the 857 selected pixels. For better visualization 
we have interpolated the data to a common wavelength grid and 
we have changed the sign in Stokes V so that the of lobe is always 
negative. All profiles are normalized to the average continuum 
intensity of the quiet Sun /qg. 

It is important to remind here that the LI 2- 2 observ- 
ing mode does not record the linear polarization profiles 
(Q, U ; see Sect. 2). In the case of strong magnetic fields 
it is plausible to recover, through the magneto-optical 
effects, the azimuth of the magnetic field vector </> from 
only Stokes / and V (Ruiz Cobo & del Toro Iniesta 
1992). However, it is unclear whether this is also the 
case when the circular polarization signals are weak 
{V/Iqs < 0.03; see Fig. 1 and also Borrero & Kobel 
2011), and the spectral resolution is limited. Therefore, 
we do not invert the azimuth angle of the magnetic 
field vector (/'(ts), and instead we fix it at a value of 
zero. Finally, we note that IMaX's spectral transmission 
profile is fully considered in the inversion. This is done 
by convolving the theoretical Stokes vector, before it is 
compared to the observed Stokes vector at each iteration 
step, with the instrument's transmission curve (see 
Fig. 2 in Martinez Fillet et al. 2011a). 

In order to improve convergence, and to reduce the 
chances of the Levenberg-Marquardt algorithm falling 
into a local minimum, each of the 857 selected pixels is 



5. INVERSION RESULTS 

The inversion of the Stokes I and V profiles from 
the 857 selected profiles provides (among other physical 
parameters) the stratification with optical depth of 
the temperature T[t^), line-of-sight velocity 1^los('''5), 
magnetic field strength B{tq), and inclination j{t^) in 
all selected pixels. Searching for similarities among all 
available results turns out to be a difficult task, as the 
inferred stratifications do not seem to follow, at first 
glance, any particular pattern. However, if one looks 
at the optical-depth dependence of the line-of-sight 
component of the magnetic field vector = Bcosj, 
one realizes that almost all pixels show a change in 
sign from B|| > to B|| < (or viceversa) at some 
optical-depth point T5 (height in the atmosphere). This 
observation allows us to classify the different results as 
a function of the optical-depth where the polarity of 
the magnetic field reverses. In particular we distinguish 
three cases: polarity change at around logrs w — 1, 
logT5 « —2, and finally, possible polarity change at 
logT5 < —3. Hereafter, these families of solutions will 
be referred to as family 1, family 2^ and family 3, 
respectively. Figures 4, 5 and 6 show the individual 
results as a function of rs from the inversion of all 
pixels belonging to each family (black-dashed lines). 
Each of these figures show: the line-of-sight component 
of the magnetic field vector (top- left panel), the 

^ In fact, as mentioned in the first paragraph in Sect. 3, this 
already happened in the V5-6 data (see Borrero et al. 2010 for 
details). 
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temperature T (top-right panel), and the hne-of-sight 
velocity Vlos (bottom- left panel). For the latter two 
physical parameters, T{t^) and 14.os(t5), we also show 
in solid-red lines the average stratification obtained 
from all the pixels belonging to a given family. For 
comparison purposes, the top-right panels in Figs. 4, 5 
and 6 also display the average temperature stratification 
(blue-solid lines) in granules (Borrero & Bellot Rubio 
2002). This model is chosen because, as mentioned in 
Sect. 3, these event occurs typically at the center or 
edges of granular cells (see also Fig. 2). In the following, 
we will describe each of the aforementioned families 
separately. 

5.1. Family 1: polarity change at logrs f« — 1. 

This family comprises 123 pixels out of the 857 
selected ones (14.3 %). As imposed by our classification 
criterion, changes sign at around logrs « —1 (Fig. 4; 
top-left panel). The temperature shows an enhancement 
of about 400 — 600 K in the mid- and upper-photosphere 
(logrs G [—1.5,-3]) with respect to a typical granule 
(Fig. 4; top-right panel). The line-of-sight velocity 
(Fig. 4; bottom-left panel) displays variations from 
extreme downflows (Vlos ~ 12 km s~^) in the upper- 
photosphere (logrs ~ —3) to large upflows in the 
mid-photosphere (Vlos ~ ^7 km s~^ at logrs ~ ^2), 
and then back to downflows in the deep-Photosphere 
(^los ~ 3 km s~^ at logrs ~ 0). Since the speed 
of sound in the solar photosphere is about Vs — 7 km 
s~^, the inferred line-of-sight velocities are close to 
supersonic. Once we consider that Vlos is only a lower 
limit of the total modulus of the velocity vector, the final 
velocities are likely to be much larger, hence supersonic. 
An additional feature is the fact that the line-of-sight 
velocity remains close to zero where the polarity of the 
magnetic field changes (logrs » — 1). 

5.2. Family 2: polarity change at logrs « —2. 

This family contains 434 of the 857 selected pixels 
(50.7 %). As imposed by the classification criterion, 
changes sign at around logrs ~ —2 (Fig. 5; top-left 
panel). Similarly to family 1, the temperature in family 
2 also shows enhancements (150 — 200 K) compared to 
an average granule. Although not as large as in the 
first case, the increase occurs over all optical depths 
(Fig. 5; top-right panel). The line-of-sight velocities are 
again large, although in this case they always involve 
upfiows (Fig. 5; bottom- left panel). These upfiows are 
visible both in the upper-photosphere (Vlos ~ —2 km 
s~^ at logrs ~ —3), and in the deep-photosphere 
(Vios ~ —7 km s~-^ at logrs ~ 0). Again, the 
line-of-sight velocity remains close to zero where the 
magnetic field changes polarity (logrs ~ —2). 

5.3. Family 3: polarity change logrs < — 3 ? 

There are 300 pixels in this family, which corresponds 
to 35.0 % of the total number. As our selection criterion 
imposes, there is no real change in the polarity of the 
magnetic field vector (Fig. 6; top-left panel). Interest- 
ingly, B|| decreases as log rs decreases. If the decreasing 



trend continues towards higher photospheric layers 
the polarity would eventually switch, although that 
would happen close to the temperature minimum. As it 
happened in the two previous families, the temperature 
in the mid-photosphere (logrs € [—1,-2]) is enhanced 
with respect to the average temperature of a granule 
(Fig. 6; top-right panel). In this case, the enhancement 
is the largest (~ 1000 K) of the three studied families. 
Finally, the line-of-sight velocity changes from large 
upflows in the mid-photosphere, Vlos ~ —8 km s~^ at 
logrs ~ ~2 (Fig. 6; bottom-left panel), to extreme 
downflows in the low-photosphere, Vios ~ 15 — 20 km 
s^^ at logrs ~ 0. As in all previous families, the 
linc;-of-siglit velocity drops to zero as the line-of-sight 
component of the magnetic field vector vanishes, which 
now happens in the upper-photosphere (Vlos — >^ at 
logrs « -3). 

6. DISCUSSION ON MODEL CHOICE 

In this section we discuss the choice of model and 
free parameters in the inversion described in Section 
4. To make this choice we look at the observed Stokes 
profiles and try to establish the most suitable model to 
fit them. The first step is to realize that, as illustrated 
in Fig. 3, all 857 selected Stokes V profiles feature the 
almost complete lack of one of the lobes (produced by 
the AM = ±1 transitions in the Zeeman pattern) in the 
circular polarization. 

The question is now to decide whether these kind of 
profiles are really asymmetric, or on the other hand, 
they are symmetric but the missing lobe in Stokes V 
is located at AA < —200 mA, hence lying just outside 
of the region scanned by IMaX. This question can be 
answered by looking at the wavelength region around 
AA » 150 - 200 ml in Figs. 4, 5 and 6 (bottom-right 
panels). The average observed Stokes profiles (filled 
circles) for all families around this region is negative 
and, although small, clearly above the noise level. We 
have established that this corresponds to the Stokes V 
signal from the nearby Fe I 5250.653 A spectral line 
by repeating all the inversions described in Section 
4, but excluding this spectral line (see Table 1). In 
this case, the fitted profiles in the bottom-right panels 
in Figs. 4-5-6 (solid-black lines) fail to reproduce the 
negative values of Stokes V around AA « 150 — 200 
mA. With this in mind we now move to AA < —200 mA 
and conclude that, if there is a missing positive lobe in 
Stokes V there, then the signal from Fe I 5250.653 A at 
AA K, 150 — 200 should also be positive. However, it is 
negative, and therefore we conclude that the observed 
Stokes V profile from Fe I 5250.653 A are indeed 
asymmetric and that there is no missing positive lobe at 
AA < -200 mA. 

Once established that the observed circular polariza- 
tion profiles are highly asymmetric, we follow Solanki & 
Montavon (1993) and Lanclolfi & Landi Degl'Innocenti 
(1996) and consider that those asymmetries are caused 
by the simultaneous effect of gradients in the line-of- 
sight velocity VloS) and in line-of-sight component 
of the magnetic field vector Bk = Bcosj. The large 
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Family 1 ; Number of pixels = 123 




Fig. 4. — Top-left panel: line-of-sight component of the magnetic field vector as a function of the optical depth By (rs). Top-right panel: 
temperature as a function of the optical depth T(r5). Bottom-left panel: line-of-sight velocity as a function of the optical depth MosCts)- 
Dashed-black lines show the results from the inversion of the Stokes profiles / and V from each of the 123 selected pixels that belong to 
family 1. The red-solid line shows the average obtained from the individual results. Blue-solid line in the top-right panel corresponds to the 
temperature stratification in the granular model by Borrero &; Bellot Rubio (2002). Bottom-right panel: average of the observed (circles) 
and of the fitted (solid line) Stokes V profiles in all pixels belonging to family 1. The mean value of from all the individual inversions 
is also indicated. 



number of nodes allowed in Vlos('''5) (see Sect. 4) is a 
direct consequence of the need to fit these extremely 
asymmetric Stokes V profiles. 

Gradients in i^y have been included by, as already 
mentioned, allowing two nodes in ^{t^) and one node in 
B{t^). This combination is more general than allowing 
two nodes for B(t^) and only one for 7(15). This is 
a consequence of the modulus of the magnetic field 
vector being defined as a positive quantity which makes, 
unlike the former case, the inversion with only one node 
in 7(''"5) unable to yield solutions where changes 
sign with optical depth. The possibility of obtaining 
solutions where i3y changes sign does not imply that all 
inferred stratification will actually have it (i.e. family 3; 
see Sect. 5.3 and Fig. 6). 

Since our discussion in the next section will rely 
heavily on the presence of a reversal in the polarity of 
the magnetic field, it is crucial to establish whether this 



feature is really needed to reproduce the observed pro- 
files. To this end we once more repeated our inversions 
(as described in Sect. 4) but employing this time one 
node in ^{t^) and two in B{t^). In this case, all retrieved 
stratification in _B||(r5) show the same sign at all optical 
depths (as in family 3). Interestingly, the quality of the 
fits worsens significantly. The average value of the merit 
function, , doubles in those pixels that belonged to 
families 1 and 2. Meanwhile, also increases in those 
pixels that belonged to family 3, but comparatively less 
(< 50 %). These results suggest that allowing two nodes 
in 7(r5) is necessary to successfully fitting the observed 
Stokes profiles, and therefore also that, the inferred 
reversal in the polarity of the magnetic field i3y (in fam- 
ilies 1 an 2) is not an artifact imposed by our choice of 
model, but rather a characteristic feature of these events. 

Finally, as previously discussed, a model that contains 
only one component assumes that the solar photophere 
is laterally homogeneous within each observed pixel. 
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Fig. 5. — Same as Figure 4 but showing the 434 pixels belonging to family 2. 



or at least, assumes that the vertical variations in the 
physical parameters play a more important role than 
the horizontal ones in the formation of the observed 
spectral line. While there is no guarantee that this is 
indeed the case, the high spatial resolution achieved by 
SuNRiSE/IMaX makes this approximation a reasonable 
first step in the study of these events. 

7. DISCUSSION AND CONCLUSIONS 

The results presented here confirm our initial con- 
clusions in our previous studies (Borrero et al. 20f0, 
2012) on these extremely shifted polarization signals. 
Namely that, a) they occur mostly at the center or 
edges of granular cells; b) they are characterized by 
supersonic upward velocities; c) they involve magnetized 
plasma, and d) magnetic fields of opposity polarities are 
oftentimes (~ 70 % of the cases) seen in their proximity 
(~ 2"). In addition to this, the inversion of the Stokes 
profiles reveals that these events seem to belong to 
three distinct families. These families frequently present 
features such: e) temperature enhancement of a few 
hundred Kelvin in the mid-photosphere; f) shift from 
supersonic upflows to supersonic downflows at some 
height in the photosphere; and g) presence of a reversal 
in the polarity of the magnetic field vector also at some 



height in the photosphere at the exact location where 
the event occurs. 

Owing to their common features, and under the 
assumption that only one physical mechanism is respon- 
sible for all the observed events, it would be almost 
straightforward to consider magnetic reconnection as 
their probable cause: magnetic field lines of opposite 
polarity coalesce and the energy stored in the magnetic 
field is released into kinetic and thermal energy. The two 
different polarities would channel the plasma in different 
directions giving rise to both positive and negative 
line-of-sight velocities (Rezaei et al. 2007; Cameron et 
al. 2011). Unfortunately, not all investigated pixels 
share the aforementioned properties. For instance, only 
some cases {families 1 and 3) show both positive and 
negative supersonic line-of-sight velocities VloSi while 
family 2 posseses only Vlos < (upflows). In addition, 
the reversal in the polarity of the magnetic field vector 
is not always present (e.g family 3; see Fig. 6), and only 
in the case of family 1 (14.3 % of the cases) the reversal 
in _B|| occurs at the same location as the change in the 
sign of ^LOS • 

On the one hand, taking into account that the inter- 
action between magnetic fields and granular convection 
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Fig. 6. — Same as Figure 4 but showing the 300 pixels belonging to family 3. 



leads to a rich variety of phenomena, it is conceivable 
that the differences between the inferred families are 
caused by the underlying physical mechanism being 
different in each case. Although there are many possible 
candidates, a search across the available literature 
(Steiner et al 1998; Cheung et al. 2008; and references 
therein) does not reveal any mechanism that reproduces 
the observational features, neither in general nor of 
the individual families, of the events studied in this 
work. For instance, the supersonic flows predicted by 
Cattaneo et al. (1990) and later observed by Ryback 
et al. (2004) and Bellot Rubio (2009), occur above 
granules and involve supersonic flows, but they are 
mostly horizontal and therefore their contribution to 
V-Los is unlikely to be large. Flux-emergence processes 
described in Cheung et al. (2008) take place also in 
granules, but they do not seem to involve very large 
upflows. Swaying motions in flux tubes (Steiner et al. 
1998) excite up-ward propagating shock fronts, but they 
occur mainly above intergranular lanes. Finally, vortex 
tubes that were originally found in granules (Steiner et 
al. 2010) have been recently associated with fast upflows 
but on nearby dark lanes (Yurchyshyn et al. 2011) and 
therefore, they probably correspond to a different kind 
of event. 



On the other hand, one could attempt to salvage 
the hypothesis of reconnection by adopting different 
views. For instance, we could argue that one cannot 
expect all families to be fully consistent with the classic 
picture of magnetic reconnection, because they might 
correspond to different stages in the temporal evolution 
of the events (see Cameron et al. 2011). In order to rule 
out or to confirm this possibility, one would need an 
uninterrupted, and possibly longer, time-series of L12-2 
data. Hopefully, this will be possible in the incoming 
second flight from SuNRiSE/IMaX that is scheduled to 
take place in the summer of 2013. It can also be argued 
that, even if events belonging to family 3 do not show a 
polarity reversal in the magnetic field, this reversal can 
indeed take place in the upper-photosphere (logrs < —3; 
see Sect. 5.3). Moreover, even if the polarity reversal is 
not present on the same pixel, in Borrero et al. (2010) 
we had already detected opposite polarities within 2" 
in 70% of the events. In the future, it would be very 
interesting to combine IMaX observations with data 
from the upcoming IRIS mission, to study a possible 
relationship between these reconnection events and 
the presence of mostly unipolar regions (coronal holes) 
and/or type II spicules in the Chromosphere (Mcintosh 
et aL 2011). 
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